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Abstract. 

We discuss two possible scenario for the direct cascade in two dimensional turbulent systems in presence of 
friction which differ by the presence or not of enstrophy dissipation in the inviscid limit. They are distinguished 
by the existence or not of a constant enstrophy transfer and by the presence of leading anomalous scaling 
in the velocity three point functions. We also point out that the velocity statistics become gaussian in the 
approximation consisting in neglecting odd order correlations in front even order ones. 
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Two-dimensional turbulence has very peculiar properties compared to 3d turbulence, see eg. ref.||l|, 
for references. As first pointed out by Kraichnan in a remarkable paper 1), these open the possibility for 
quite different scenario for the behaviors of turbulent flows in 2d and 3d: if energy and enstrophy density 
are injected at a scale Li, with respective rate e and ~ 'SL~^, the 2d turbulent systems should react such 
that the energy flows toward the large scales and the enstrophy towards the small scales. One usually refers 
to the infrared energy flow as the inverse cascade and to the ultraviolet enstrophy flow as the direct cascade. 
In the (IR) inverse cascade, scaling arguments lead to Kolmogorov's spectrum, with E{k) ~ e^/^ k~^^^ for 
the energy and (Su) ~ (er)^/'^ for the variation of the velocity on scale r. In the (UV) direct cascade, scaling 
arguments give Kraichnan's spectrum with E(k) ^ e^{^ k^^ for the energy and {Su) ^ (cmi?'^)^^'^ for the 
velocity variation. The aim of this Letter is to discuss a few possible scenario for the influence of friction on 
the direct enstrophy cascade of 2d turbulence. The picture which emerges is that at small scales the velocity 
field has a smooth gaussian component, with amplitude diverging in the frictionless limit, supplemented by 
possible anomalous corrections whose characters depend whether there is entrophy dissipation or not in the 
inviscid limit. We also point out that friction provides a way to regularize amplitudes in the inviscid limit. 

A few experimental facts. Two-dimensional turbulence has recently been observed in remarkable ex- 
periments providing new informations on these peculiar characteristics. See refs.[^ ^ for a precise description 
of the experimental set-up. The turbulent flow takes place in a square cell such that the injection length is 
Li ~ 10 cm. The bottom of the cell induces friction with a coefficient r ~ 25 s. The dissipation length 1,1 is 
of order 1 mm and the UV friction length If of order 0.5 cm. (r, l^, If are defined below). The enstrophy 
transfer rate determined from the three point structure function has been evaluated in |^ as 77^ ~ 0.4 s~^. 
The direct cascade is observed at intermediate scales between 1 cm and 10 cm. Besides the absence of any 
experimentally significant deviations from Kraichnan's scaling, two remarkable facts have been observed: (i) 
the vorticity structure functions are almost constant on the inertial range with the two-point function of 
order 10 s^^; and (ii) the vorticity probability distribution function is almost symmetric and almost, but 
not quite, gaussian. 

A model system. As usual, to statistically model turbulent flows we consider the Navier-Stokes equation 
with a forcing term. Let {x, t) be the velocity field for an incompressible fluid, V-u = 0. The Navier-Stockes 
equation with friction reads: 

dtu^ -t- {u ■ - uV^u^ + -u^ = -V^p + f' (1) 

T 

with p the pressure and /(x, t) the external force such that V • / = 0. The pressure is linked to the velocity 
by the incompressibility condition: V^p = — V^V-'u'm-'. The friction term ^u^ is introduced in order to 
mimick that in physical systems the infrared energy cascade terminates at the largest possible scale. The 
vorticity w, with w = CijdiUj, is transported by the fluid and satisfies: dtu> + {u ■ V)u; — vSl'^io + uj /t ^ F 
with F ~ £ijdifj 1^. We choose the force to be gaussian, white-noice in time, with zero mean and two point 
function: {f^{x,t) f''{y,s)) = C^'^{x — y) 5{t — s) where C^'^{x), with C^^[x) = 0, is a smooth function 
varying on a scale Li and fastly decreasing at infinity. The correlation function of the vorticity forcing term 
is: {F{x, t) F{y, s)) = G[x — y) 5{t — s) with G — — V^C. We shall assume translation, rotation and parity 
invariance and expand C^^{x) as: C^^{x) = eS^^ -Iw (3r^i5*-' - 2x^x^)/8 + ■ ■ ■ with — x^Xk- The scale Li 
represents the injection length. The inverse cascade takes place at distances Li <^ x <^ L f ~ r'^/^e^/^. At 
finite viscocity, there are two ultraviolet characteristic lengths, the usual dissipative length Id — v^^'^lw^^^ 
and another friction length If ~ v^/"^ r^/^ above which friction dominates over dissipation. The direct cascade 
takes place at scale If <^ x Li. We shall consider the inviscid limit at fixed friction, ie. — > at r fixed. 

Note that (/// Idf — TeU^ is large in the quoted experiments. 

The fundamental property of two dimensional turbulence recognized by Kraichnan |^ and Batchelor ||^ is 
that the energy cascades towards the large scales because it cannot be dissipated at small scales. In absence 
of friction the system does not reach a stationnary state. But friction provides a way for the energy to escape 
and for the system to reach a stationnary state. The fact that the energy is not dissipated at small scales 
translates into the vanishing of the averaged energy dissipation rate, ie. i/((Vw)^} = in the inviscid limit, 
which means that there are no energy dissipative anomaly |8|, ie. 

limi.(Vu)^,) =0 (2) 
inside any correlation functions. However there could be enstrophy dissipative anomalies in the sense that 
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i^(Vaj)^ does no vanish in the inviscid correlation functions: 



e^(x) = lim i/(Vw)L) (3) 

Below we shall discuss possible consequences of the vanishing or not of in presence of friction. 

Two and three point velocity correlation functions. In absence of energy dissipative anomaly, 
the mean energy density relaxes in the inviscid limit according to dt{\) + ^{y?) = e, showing that e is 
effectively the energy injection rate. It reaches a stationary limit with {v^) = er. Similarly the enstrophy 
density evolves according to dt{ijj'^) + 2i^((Vw)^) + -(w^) — 2ew showing that Iw is the enstrophy injection 
rate. Let = lim^^o ^ii'^^)'^) be the enstrophy dissipation rate, e^, < Zw In the stationary limit at finite 
friction the enstrophy density is finite and equal to: 

= r (e^ - e„) (4) 

le. the enstrophy density is equal to the difference of the enstrophy injection and the enstrophy dissipation 
rates times the friction relaxation time. The vorticity two point correlation function {uj{x)u!{0)) stay finite 
since it is bounded by (w^). Hence the vorticity correlation cannot diverge as a: ^ 0, and it cannot have 
a negative anomalous dimension. Scalings {u!{x)uj{Q)) ~ r~'> or (logr)^ with ^' > are forbidden in 
presence of friction. At short distance one then may have: 

{lu{x)lo{0)),^o =n- {Te.^)A {r/Lf^^ + ■ • ■ (5) 

with ^2 > 0. We shall denote the limiting value hy Q. = T^e^ — Vw)- By eq(^) below, rj^ will be identified 
as the enstrophy transfer rate. The amplitude A, rj^ and as well as the anomalous exponant ^2, are 
functions of the dimensionless parameter r'^e^ . 

Since V^(((5m)^) = 2{uj(x)ijj{Q)), finiteness of the vorticity two point function at coincidant points imply, 
{{Suf')) ~ ^r^, but the subleading terms could have anomalous scalings: 

— . _2 /3 

The leading scaling is the normal Kraichnan scaling but the amplitude is different, ie. fl instead of . There 
is a crossover at scale rc, with r^^^ ^ (il/re^), above which the second anomalous contribution dominates. It 
is interesting to compare the occurence of both the smooth term in eq. and the anomalous contribution 
r^^^ in eq.(^ with the passive scalar problem with friction |^ in which anomalous zero modes occur in 
the passive scalar correlations only if the velocity flow is regular. The smooth term in {(Su)"^) does not 
contribute to the asymptotic large k behavior of the energy spectrum which thus scales as fc~^~^^^ . Namely: 

E{k) ^ -^!^!^;i|±iH) (re^)Ak^' (kL)-'^^ (7) 

as — > 00. The exponant ^2 is expected to be universal but not the amplitude A. 

As usual, stationarity of the two point correlation functions gives in the inviscid limit, for x ^ 0: 

V';{{Su''){Suf) + 2{{Suf)/T = 2(2e-C(x)) ~e^r2 + ... (8) 

It allows us to determine exactly the three point velocity inviscid structure functions: 

3Ae,,, 



4(6 + i)'(6 + 2)(6 + 3) 
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Without dissipative anomaly at finite friction. If there is no enstrophy dissipative anomalies in 
presence of friction, 'iw — and hence 77^, = assuming continuity of the correlations. The three point 
velocity structure (bI) has an anomalous scaling at leading order. It behaves as 



r 
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since the first term rj^r^ vanishes. As a consequence the enstrophy transfer is not constant through the 
scales, ie. there is no enstrophy cascade. The anomalous scaling manifests itself only at a subleading order 
in the velocity two point function but the enstrophy structure function is anomalous: 

iiScuf) - 2(Te^)A(r/i)2«^ 

Moreover in absence of enstrophy dissipative anomaly the one point vorticity correlations are gaussian, 

(exp(sw)} = exp(s2(re^)/2) (10) 

This is a direct consequence of the stationarity of the vorticity correlations. It is not true for the higher point 
functions. It could provide a way to check whether the enstrophy dissipation vanishes or not in presence 
of friction. Below we shall present a perturbative argument in favor of the absence of enstrophy dissipative 
anomalies in presence of friction. 

In this scenario with t^, = 0, the quasi symmetry of the probability distribution may be explained by the 
fact that the even order correlations are parametrically larger that the odd order correlations, as it is for 
the two and three point function, cf eqs.(^,^ with rjw = 0. Namely the even order correlation have leading 
normal but anomalous subleading contributions whereas the odd order correlations would have anomalous 
contributions at leading order which decrease faster at small scales. 

A scenario for the frictionless limit. Experimentally, (r'^e^) is large and we are thus interested in 
the frictionless limit. We shall discuss this limit assuming the absence of dissipative anomaly at r finite. 
Demanding that the three point structure function reproduces the known frictionless exact result with 
{{6u)^) ~ r'^ as r oo imposes that the anomalous exponant ^2 vanishes and that the amplitude A goes to 
one as T — > CX3. Moreover, demanding that the velocity structure function is finite as r — + 00 requires: 

6 = C2(T3e»)-'/'+o((T3e^„)"'/') and A = 1 + c{T'e^)-'/' + o{{T'e^,)-'^'') (11) 

with C2 and c constant. Since ^2 is supposed to be universal but not the amplitude A, we expect ^2 universal 
but not c. As a consequence, at fixed non zero distance r, the vorticity correlation behaves logarithmically 
as: 

{u;{xMO)) ^el/^ [c-2C2log(r/L)]+o((r3e^)-V3) (12) 

Here we assume that the terms represented by the dots in eq.(^) which are subleading, and thus negligeable, 
at finite friction do not become predominant in the frictionless limit. But this is the simplest scenario. 

The vorticity structure function, estimated as {{Suj)'^) — ItZw—'^w'^ (c— 2C2 log(r/L)), is then dominated 
by the constant term since the enstrophy density diverges as r — s- 00. This may explain the experimental 
fact that vorticity correlations are almost constant in the inertial range. This scenario is compatible with 
experimental data since experimentally (((^tj)^) ~ 10 s^^ whereas 2Te^ ~ 20 with an estimated error of 
20*^/0 to 30°/o. The velocity structure function will then be: 

{{8uf) ^ 6^/3 [c/2 + C2 - C2 log(r/L) + ■ • •] 

The smooth is not universal but the logarithmic correction is (since C2 is expected to be universal). In 
the same way the energy spectrum (Q) Eik) becomes proportional to k^^ with: 

i?(fc):^2C2e^/'fc-' (13) 

with Kraichnain's constant Ck = 2^2 which is universal if the anomalous dimensions at r finite are universal. 
Note that the smooth but non universal terms in the velocity correlation does not contribute to the energy 
spectrum at large momenta. Of course, there could be logarithmic correction to the spectrum if ^2 — 0. 

With dissipative anomaly at finite friction. If the enstrophy dissipation is non zero in presence 
of friction, the enstrophy density is finite but smaller than the injected enstrophy density: (w^) < re^,. 
The three point velocity structure function then scales as (((5m)'^) ~ r/^r'^ which allows us to identify 77^ as 
the enstrophy transfer rate, ie. there is enstrophy cascade with transfer rate r\w The quasi symmetry of 
the probability distribution function would hold only if the entrophy transfer rate is much smaller that the 
entrophy injection rate r/™ <C eu,. There are many scenario for the frictionless limit depending on how riw/f-w, 
which is a function of r^e^, behaves as r — > 00. 
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Neglecting odd order correlations. We now discuss what would be the consequences of neglecting 
the odd order correlations in front of the even order ones. This is suggested by the quasi symmetry of the 
experimentally measured probability distribution. Consider the equations encoding the stationarity of an 
even number of velocity differences, ie. dt{{Su^^) ■ ■ ■ {6u^^")) = 0. Both the terms arising from the advection 
{u ■ V)m or from the pressure Vp involve correlation functions with an odd number of velocity insertions. 
Assuming that the odd order velocity correlations are much smaller that the even order ones, these terms 
are much smaller and negligeable compared to the remaining terms which are those arising from the friction 
and from the forcing and which involve an even number of velocity insertions. Hence at the leading order 
the stationarity condition of the even order correlations reads: 

. . - ^ 

— {{Su'')- ■ ■ ((5m*^"))("^ = 2^2(6^11]'') {■ ■ ■ i^u^') ■ ■ ■ i^u'") ■ ■ (14) 

p<q 

with SC'^^ (x) — C*-' (0) — C'^^x). The overhatted quantities are omitted in the correlation functions. In 



eqs.(|14[) we have used the fact that there is no energy dissipative anomalies. The solution of eqs.(14) is 
provided by a gaussian statistics with zero mean and two-point function ~ ^^^^^{Sr^S''^ —2x^x^ ). 

Hence, in this approximation the velocity structure functions scale as: 

F^°J = ((5u)2"}(0) ~ const, (re^r^)" (15) 

This is compatible with eq.(|^) since fl — re^ if 77^, To be consistent it also has to solve the stationarity 

equations for the odd order correlations which at this order read: 

with V^p = —S/^V^u^u^ . Since the statistics is gaussian it is enough to check them for n — 2. The gradiant 
of the pressure has scaling dimension one. Then the l.h.s. up to O(r^) is a rank three symmetric 0{2) 
polynomial tensor, transverse along all coordinates, invariant by translation and of scaling dimension three. 
There is no such tensor and thus the l.h.s. is 0{r^) for n = 2. In other words, the n-point velocity correlations 
computed using the gaussian statistics are zero modes. 

In view of the exact formula eqs.(^,^, this approximation would be better for the transverse velocity 
correlations. Note also that using this approximation to compute the one point vorticity functions (w") gives 
the exact result ( p^ ) if there is no enstrophy dissipative anomaly. Of course there would be deviations to 
this approximation: 

with corrections (5F„ including both contributions with normal scaling or with possible anomalous scaling as 
in eq.(|^,^. In the scenario without enstrophy dissipation, these anomalous contributions will be subleading 
in the even order correlation but dominant at small scales in the odd order correlations. On contrary if 
the enstrophy transfer is not vanishing these corrections include contributions with normal scaling, since 
5F2 — rji^rr"^ and SF3 ~ rj^r^ as follows from eqs.(^,^. Stationnarity of the three point functions leads 
then to an equation for the four point correlations of the form \/'^F4 + -F3 — with the transverse 
gradiant. Using ^3 ~ rj^r'^ and naive scaling, ie. omitting possible zero mode contributions, gives a correction 
6F4 ~ 2^7-4 ^]^g gj.g^ order gaussian contribution i^j"' which satisfies V^F^^^ = 0. Similarly, stationarity 
of the four point functions now give equations of the form V^Fs + ^SF^ ~ rj^Te^r^. For t^Iw ^ 1 as in 
one may neglect ■^SF4 in front of rj^iTlwr^ , so that 5F^ ~ rjwTlwT^ ■ Recursively, naive dimensional analysis 
apply to stationarity equations with t^Iw ^ 1 yield to first order in 77^,: 5Fn+2 — Jytu t^~" (T'^e^,)'^ r"+^ 
with k the integer part of n/S. Of course 5F2n ^ are small for r^tw 3> 1 > {rjw/^w)- But zero modes 
of the stationarity equations which have been neglected in this dimensional analysis may give additional 
contributions to 5Fn. These zero modes will be the leading terms for the odd order correlations if the 
enstrophy dissipation and thus the enstrophy transfer rate vanish. 

MSR formalism. We now discuss how these scenario could fit in a field theory approach to the 
problem. The corrections ( pT|) of the anomalous exponant as a function of the friction coefficient have a 
natural interpretation in the MSR formalism |ll0| . The MSR formalism provides a way to compute the 
inviscid vorticity correlations using a path integral with action 

S = i 1 dtdxip(^^^t){dti^ + u-Vuj + -u;)(^^^t) + 7: [ dtdxdy ip(^^t)G(a:-y)V(y,t) (17) 
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Note that this action describes the inviscid hmit v = 0. 

Let us first define Lu{x,t) = ew'^Lu{teU'^ , x) and (p{x,t) — e-ui^^'^ipiteU^ ,x). This corresponds to select 

_— 1/3 

configurations for which the dominant time scales and amplitudes are ■ It allows us to present the 
theory as a perturbation of the frictionlcss action Sq: 

S = So+^9-'^-'Jdtdx^,^,,^,^,, (18) 

with g = T^Iw as dimensionless expansion parameter. This shows that the pertubation around the frictionless 
theory is controled by (r^e^)"^/'^ as in eg. (p^ . Of course implementing this perturbation remains a challenge 
as the unperturbed theory is unkown. 

Alternatively let us define Lj{x,t) = {tZwY^'^^^' {t/r^x) ^i^d (p{x,t) = {Tliu)^^^^tp'{t/T,x). This corre- 
sponds to select configurations with time scale r. The theory is then formulated as a perturbation of the 
gaussian action by the advection term with g^/^ = (r'^e^)^/^ as coupling constant: 

S = 5'gauss + ig^^^ J dtdx ^) 4) • V)W(^ (19) 

where S'gauss is the gaussian action obtained by neglecting the advection term. This is the statistics we 
found by neglecting odd order correlation functions. In the unperturbed theory (w^) = re^, and there is no 
anomaly. Perturbation theory is then computable. We checked that there is no occurence of the dissipative 
anomaly in first order (one loop) in the perturbation theory. This is linked to the fact that the friction 
regularizes the theory in the inviscid limit such that there is no divergence in the Feymann amplitudes, at 
least in the first order we checked. This is expected to be true to all orders. In other words, there is no 
pertubative contribution to the enstrophy transfer rate, ie eu, = 77^, — in perturbation theory. 

Non perturbative contributions may arise due to instanton contributions to the path integral. Let us for 
exemple define uj{x, t) — T~^uj{t/T, x) and (p{x, t) — {T^e.u])~^^{t/T, x). It allows us to write the action as S* = 
{T^f^ui)~^S with S dimensionless. In the saddle point approximation, instantons could then potentially give 
a non perturbative contribution to the enstrophy transfer rate with 77^ ~ (const. )t~^ exp(— const. /(t'^Eu,)). 
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